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ABSTRACT 

In this study we investigate the relationship between the star formation rate (SFR) 
and AGN luminosity (Lagn) for ~2000 X-ray detected AGN. The AGN span over 
three orders of magnitude in X-ray luminosity (10^^ <L 2 _ 8 keV< 10^®’^ erg s“^) and 
are in the redshift range 2 ; = 0.2 - 2.5. Using infrared (IR) photometry (8 - 500/im), 
including deblended Spitzer and Herschel images and taking into account photometric 
upper limits, we decompose the IR spectral energy distributions into AGN and star 
formation components. Using the IR luminosities due to star formation, we investigate 
the average SFRs as a function of redshift and AGN luminosity. In agreement with 
previous studies, we find a strong evolution of the average SFR with redshift, tracking 
the observed evolution of the overall star forming galaxy population. However, we find 
that the relationship between the average SFR and AGN luminosity is broadly flat 
at all redshifts and across all the AGN luminosities investigated; in comparison to 
previous studies, we find less scatter amongst the average SFRs across the wide range 
of AGN luminosities investigated. By comparing to empirical models, we argue that 
the observed flat relationship is due to short timescale variations in AGN luminosity, 
driven by changes in the mass accretion rate, which wash out any underlying corre¬ 
lations between SFR and Lagn- Furthermore, we show that the exact form of the 
predicted relationship between SFR and AGN luminosity (and it’s normalisation) is 
highly sensitive to the assumed intrinsic Eddington ratio distribution. 

Key words: galaxies: active; — galaxies: evolution; — X-rays: galaxies; — infrared: 
galaxies 


1 INTRODUCTION 

One of the key outstanding problems in studies of galaxy 
evolution is understanding the connection between active 
galactic nuclei (AGN) and star formation. Both AGN activ¬ 
ity and star formation are predominately dependent on the 
availability of a cold gas supply from the galaxy, as it is the 
fuel of both processes, and therefore a first order connection 
between these two processes may be expected. However the 
scales of AGN activity and star formation are very different, 
which has lead to suggestions that any tight connection be¬ 
tween them must be due to one process regulating the other 
(see Alexander & Hickox 2012, Fabian 2012, and Kormendy 
& Ho 2013 for recent reviews). 


There are several pieces of empirical evidence for at 
least a broad connection between AGN activity and star 
formation. For example, the tight correlation observed be¬ 
tween the mass of the super-massive black hole (SMBH) 
and the galaxy spheroid for galaxies in the local universe 
(e.g., Kormendy & Richstone 1995, Magorrian et al. 1998), 
serves as archaeological evidence of a connection between 
the growth of the SMBH (through mass accretion, where 
it becomes visible as AGN activity), and the growth of the 
galaxy (through star formation). Additionally, observations 
of AGN have found that the volume average of the SMBH 
mass accretion rate tracks that of the star formation rate 
(SFR), within ~3-4 orders of magnitude, up to redshifts of 
2 ~2 (e.g., Heckman et al. 2004; Merloni, Rudnick & Di 
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Matteo 2004; Aird et al. 2010) suggesting a co-evolution of 
AGN and star forming activity. Despite how significant these 
results may appear, they only provide indirect evidence for 
a relationship between AGN activity and star formation and 
cannot place strong constraints on the form of the relation¬ 
ship. 

To acquire more direct evidence on the form of the rela¬ 
tionship between AGN activity and star formation requires 
sensitive measurements of the AGN and star forming lumi¬ 
nosities of individual galaxies. X-ray and far-infrared (FIR; 
A = 30 — 500 fim) observations are ideal for quantifying the 
amount of AGN and star formation activity, respectively. 
A key advantage of X-ray observations, specifically in the 
hard band (e.g., 2-8 keV), over other tracers of AGN 
activity, is that they are not greatly affected by the pres¬ 
ence of obscuration and contamination effects from the host 
galaxy (see sections 1 and 2 of Brandt & Alexander 2015 
for more details of the use of the X-ray as an AGN tracer). 
A key advantage of FIR observations, as a measurement of 
star formation, is that they trace the peak of the obscured 
emission from star forming regions surrounded by cold gas 
and dust. Even though the FIR provides an indirect tracer 
of star formation, a significant advantage over more direct 
tracers, such as the UV and optical emission from the young 
massive stars, is that it does not suffer significantly from ob¬ 
scuration (e.g. Kennicutt 1998; Calzetti et al. 2007; Calzetti 
et al. 2010; see also section 2.2 in Lutz 2014). Indeed, as 
shown by Dominguez Sanchez et al. (2014), for luminous 
infrared galaxies (FIR luminosities of Lfir -J lO"^"^ erg s“^) 
more than 75% of the total emission due to star formation 
is produced at FIR wavelengths, a fraction that increases at 
higher Lfir.^ However, the AGN can also contribute to the 
FIR luminosity due to the thermal re-radiation of obscur¬ 
ing dust from the surrounding torus (e.g. Antonucci 1993). 
Hence, for the most reliable measurements of the star for¬ 
mation it is important to apply decomposition methods of 
the AGN and star formation components at infrared wave¬ 
lengths (e.g. Netzer et al. 2007; Mullaney et al. 2011; Del 
Moro et al. 2013; Delvecchio et al. 2014). 

A number of studies have used X-ray and FIR obser¬ 
vations to understand the connection between distant AGN 
activity and star formation by measuring the mean SFRs 
of AGN and star forming galaxy samples (e.g.. Lutz et al. 
2010; Shao et al. 2010; Mainieri et al. 2011; Mullaney et al. 
2012b; Rovilos et al. 2012; Santini et al. 2012; Harrison et al. 
2012; Rosario et al. 2013a,b; Lanzuisi et al. 2015). The main 
results shown by these studies are that: (1) the average star 
formation rates ((SFR)) of AGN track the increase with red- 
shift found for the overall star forming galaxy population; 
(2) the (SFR) of AGN are higher than those of the overall 
galaxy population (i.e., when including quiescent galaxies); 
and (3) the specific SFRs (i.e., the ratio of SFR over stellar 
mass, which serves as a measure of the relative growthrate 
of the galaxy) of AGN are in quantitative agreement with 
those of star forming galaxies. The majority of the current 

^ We note that for less luminous infrared galaxies (Lpur < 
10"^^ erg s“^) Dominguez Sanchez et al. (2014) find that the FIR 
emission accounts for ~50% of the total emission due to star 
formation. However in this work we find that our galaxies have 
average Lfir lO"^"^ erg s“^ and so the majority of the star for¬ 
mation is expected to be produced at FIR wavelengths. 


studies also find no correlation between the AGN luminosity 
and (SFR) for moderate luminosity AGN (X-ray luminosi¬ 
ties of L 2 - 8 keV< 10'*'* ergs“^; e.g.. Lutz et al. 2010; Shao 
et al. 2010; Mullaney et al. 2012b; Rovilos et al. 2012; Har¬ 
rison et al. 2012). However, there are significant disagree¬ 
ments in the results for high luminosity AGN (L 2 - 8 kev)^ 
10^"^ erg s“^). There are studies arguing that the (SFR) in¬ 
creases at high AGN luminosities (e.g.. Lutz et al. (2010); 
Rovilos et al. 2012; Santini et al. 2012), a result that seems 
in agreement with the concept of AGN and star formation 
activity being connected due to their mutual dependence on 
the cold gas supply in the galaxy. Other studies have argued 
that the SFR decreases at high AGN luminosities (e.g., Page 
et al. 2012; Barger et al. 2014), potentially suggesting that 
the AGN may be responsible for reducing or even quenching 
the ongoing star formation (a result also inferred by some 
simulations of galaxy evolution; e.g., Di Matteo, Springel 
& Hernquist 2005; Hopkins et al. 2005; Debuhr, Quataert 
& Ma 2012). There are also studies arguing that (SFR) re¬ 
mains constant up to high AGN luminosities (i.e., a broadly 
flat relationship; e.g., Harrison et al. 2012; Rosario et al. 
2012; Azadi et al. 2015), extending the trend seen for mod¬ 
erate luminosity AGN. Nevertheless, the difference in the 
conclusions of such studies could be attributed to the low 
source statistics for high luminosity AGN, and strong field 
to field variations (e.g., Harrison et al. 2012). For example, 
Harrison et al. (2012) demonstrated that when using a large 
high luminosity AGN sample the broadly flat relationship 
between (SFR) and AGN luminosity found for moderate lu¬ 
minosity AGN continues to high luminosities, with no clear 
evidence for either a positive or negative correlation (see also 
Harrison 2014 for a recent review). 

To first order a flat relationship between (SFR) and 
AGN luminosity can seem surprising, since it appears to 
suggest the lack of a connection between AGN activity and 
star formation. However, Hickox et al. (2014) have shown 
that a true underlying correlation between AGN luminosity 
and (SFR) can be masked if the AGN varies significantly 
(i.e., by more than an order of magnitude) on much shorter 
timescales than the star formation across the galaxy. In fact, 
observational studies such as Rafferty et al. (2011), Mullaney 
et al. (2012a), Ghen et al. (2013), Delvecchio et al. (2014), 
and Rodighiero et al. (2015) have shown that when the av¬ 
erage AGN luminosity is calculated as a function of SFR 
(i.e., taking the average of the more variable quantity as a 
function of the more stable quantity) a positive relationship 
is found, suggesting that AGN activity and star formation 
are correlated on long timescales. Studies using small scale 
hydrodynamical simulations of SMBH growth (e.g., Gabor 
& Bournaud 2013; Volonteri et al. 2015) have indeed sug¬ 
gested that AGN activity can vary by a typical factor of 
~ 100 over ~Myr timescales, which results in a flat rela¬ 
tionship between (SFR) and AGN luminosity over a wide 
range of AGN luminosity. These studies therefore demon¬ 
strate that the relationship between AGN luminosity and 
(SFR) can potentially place constraints on the variability 
of mass accretion onto the SMBH in galaxies. However, to 
date, the observational constraints of the (SFR) of AGN as a 
function of AGN luminosity and redshift have lacked the ac¬ 
curacy to be able to distinguish between the different SMBH 
mass accretion models. 

Most of the current studies on the (SFR) of distant 
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X-ray AGN suffer from a variety of limitations, which af¬ 
fect the accuracy of (SFR) measurements, such as: (1) small 
number of sources, which can lead to large statistical un¬ 
certainties, particularly at high AGN luminosities; (2) high 
levels of source confusion at FIR wavelengths, which can 
cause the overestimation of the flux; (3) use of a single FIR 
band from which to derive SFRs, which will result in large 
uncertainties on the (SFR) and will not take into account 
possible contamination of the SFR measurements from the 
AGN; (4) neglect of the information that can be obtained 
from the photometric upper limits of the FIR undetected 
AGN, which make up the majority of the distant AGN in 
X-ray samples (this final point is not applicable for studies 
that use stacking analyses). 

In this work we aim to overcome the limitations out¬ 
lined above by exploiting a large sample of X-ray detected 
AGN with deep and extensive multi-wavelength data, for 
which we perform spectral energy distribution (SED) fit¬ 
ting on a source by source basis, and measure the SFR for 
each source in our sample. We use deblended FIR photom¬ 
etry from Herschel, which provides the best constraints on 
the FIR fluxes of individual sources by reducing the con¬ 
tamination due to blended and confused sources, the most 
significant drawback of the Herschel field maps. Further¬ 
more, we make use of the photometric upper limits in the 
fitting procedure to better constrain the SED templates and 
SERs. We finally calculate the (SFR) values as a function of 
X-ray luminosity, with the inclusion of sources with only up¬ 
per limit constraints using survival analysis techniques (e.g., 
Feigelson & Nelson 1985, Schmitt 1985). Our methods en¬ 
sure the use of all available data (i.e. photometric detections 
and upper limits, SFR measurements and upper limits) to 
provide improved (SFR) as a function of X-ray luminosity 
and redshift. In Section 2 we outline the photometric cata¬ 
logues used in this work, as well as the choice of redshift and 
the choice of matching radii between photometric positions. 
In Section 3 we analyse our methods of SED-fitting as well 
as the calculation of the average IR luminosity due to star 
formation ((Lir^sf)). Finally, in Section 4 we present and 
discuss our results. In our analysis we use Ho = 71km s“^, 
Hm = 0.27, Ha = 0.73 and assume a Ghabrier (2003) initial 
mass function (IMF). 


2 AGN SAMPLE, IR PHOTOMETRY AND 
REDSHIFTS 

In this work we use the available Mid-IR (MIR; A « 3 — 30) 
to Far-IR (FIR; A « 30 — 500 fim) photometric data to 
constrain the average SFRs of a large sample of X-ray de¬ 
tected AGN over the redshift range « « 0.2 - 2.5. To con¬ 
struct a large sample of X-ray detected AGN we make use of 
three fields with deep X-ray observations: (1) Chandra Deep 
Field North (GDF-N; Alexander et al. 2003b), (2) Chan¬ 
dra Deep Field South (CDF-S; Xue et al. 2011), and (3) 
a combination of G/iandro-COSMOS (C-GOSMOS; Elvis, 
Civano & Vignali et al. 2009) and AMM-GOSMOS (Gappel- 
luti et al. 2009). To construct our final AGN sample we ob¬ 
tain the MIR and FIR photometry from observations of the 
X-ray deep fields made with the Spitzer (Werner et al. 2004) 
and Herschel (Pilbratt et al. 2010) space observatories. 
The recent Hersehel observational programs PEP/GOODS- 


H (Lutz et al. 2011; Elbaz et al. 2011) and HerMES (Oliver 
et al. 2012) in the three fields of GOODS-N, GOODS-S, and 
COSMOS, covering the wavelength range of 70 - 500/im are 
our main source of the FIR photometry (details in §2.2). We 
therefore restrict the CDF-N and CDF-S X-ray catalogues 
to these regions with sensitive MIR-FIR coverage, i.e. the 
GOODS-N and GOODS-S with areas of 187 arcmin^ each, 
but use the full 2 deg^ of COSMOS. In total these areas cover 
3609 X-ray sources. Figure 1 shows the X-ray sources in 
GOODS-N, GOODS-S, C-COSMOS, and AMM-COSMOS 
in the L 2 - 8 keV - « plane. 

In the following subsections we describe our sample se¬ 
lection and the catalogues used for the sample. In §2.1 we 
present the X-ray observations used to define our AGN sam¬ 
ple and to determine their X-ray luminosities. In §2.2 we 
present the MIR and FIR photometric catalogues used to 
constrain the SFRs of the AGN hosts via SED fitting. In 
§2.3 we describe the method of matching the X-ray sources 
to the MIR and FIR catalogues and the redshift counter¬ 
parts. 


2.1 X-ray Data 

To select the sample of AGN for our study we use the pub¬ 
licly available X-ray catalogues for the CDF-N (Alexander 
et al. 2003b), CDF-S (Xue et al. 2010) and COSMOS (Elvis, 
Civano & Vignali et al. 2009; Cappelluti et al. 2009) fields, 
restricted to the areas covered by PEP/GOODS-if and Her¬ 
MES observations as described above (see Figure 1). For the 
COSMOS field we use the C-COSMOS X-ray catalogue as 
the primary sample, while for the sources over the larger re¬ 
gion, not covered by Chandra, we use the AMM-COSMOS 
catalogue. Rest-frame, hard band 2-8 keV luminosities were 
calculated following Alexander et al. (2003a) with the equa¬ 
tion: 

L2-8keV= 47r X Dr X F2-8keV X (1 -f z)^'" (1) 

where F 2 - 8 keV is the observed X-ray hard band flux (2-8 
keV), Dl is the luminosity distance, « is the redshift (see Sec¬ 
tion 2.3), and F is the photon index used for k-corrections, 
which was fixed to a standard value of F = 1.9 (e.g., Nandra 
& Pounds 1994). Although the hard band observations in 
CDF-N and CDF-S are in the energy range of 2-8 keV, the 
C-COSMOS and AMM-COSMOS catalogues report hard 
band fluxes of the energy range of 2-10 keV. To convert the 
2-10 keV to 2-8 keV fluxes we assume F = 1.9 which yields 
a conversion factor of 0.85. 

For the 20% of X-ray sources in our final sample (see 
below) not detected in the hard band we used the full band 
of 0.5 - 8 keV (or the soft band of 0.5 - 2 keV if undetected in 
the full band) to estimate the hard band flux. We estimated 
the hard-band flux assuming a F = 1.4 spectral slope, unless 
this provided a measurement greater than the hard-band 
upper limit, in which case we assumed a F =2.3 spectral 
slope; the assumed range in spectral slope is motivated by 
the range observed in AGN (e.g., Nandra & Pounds 1994; 
George et al. 2000). Overall, with this procedure, the hard- 
band fluxes were estimated assuming F = 1.4 for 19% and 
r = 2.3 for 1% of sources in our sample (see Figure 1). 
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2.2 Mid-IR &L Far-IR Data 

To measure the SFRs of our AGN sample we need to reli¬ 
ably constrain the IR Inminosity due to star formation and 
remove any contribution from the AGN. To do this we need 
data covering both the MIR and FIR wavelengths for each 
source in our sample (e.g., Mullaney et al. 2011). We ex¬ 
ploit available photometry in the wavelength range of 8/rm 
- SOO^m, provided by observations carried out by: Spitzer- 
IRAC at 8/rm; Spitzer-IRS at 16/rm; Spitzer-MlPS at 24/rm, 
70/rm; Herschel-PACS at 70, 100, IGO^m; and Herschel- 
SPIRE at 250, 350, 500/rm. One of the advantages of our 
study over several previous studies, is the use of catalogues 
of deblended FIR Herschel images (details below). The de¬ 
blending of sources in the PACS and SPIRE observations 
allows us to overcome the blending and confusion issues en¬ 
countered in dense fields that can lead to an overestima¬ 
tion of the flux densities (e.g., Oliver et al. 2012; Magnelli 
et al. 2013). It also ensures the direct association between 
the measured FIR flux densities and the sources used as pri¬ 
ors in the deblending process. In addition to this, we also 
make sure that we have a reliable photometric upper limit 
for sources not detected in the FIR. This enables us to con¬ 
strain the star forming galaxy templates and gain an upper 
limit on the IR luminosity due to star formation, as we de¬ 
scribe in §3.1. 

The MIPS 24/rm photometric catalogues that we use 
were created by Magnelli et al. (2013). These catalogues 
are made by simultaneous PSF fitting to the prior positions 
of 3.6/im sources. The catalogues were limited to a dn de¬ 
tection limit at 24/im going down to 20^Jy in GOODS-N 
and GOODS-S, and 50/iJy in COSMOS. The PACS 70/rm, 
lOO/rm and 160/rm catalogues were also created by Mag¬ 
nelli et al. (2013) using the MIPS 24/im detected sources, 
described above, as the priors for the deblending of the 
PACS maps. Only sources with at least a 3a detection at 
MIPS 24/im were used as priors and the resulting PACS 
catalogues were also limited to a 3a detection limit. ^ The 
SPIRE 250pim, 350/rm, and 500/rm catalogues were created 
following the method described in Swinbank et al. (2014), 
again using these MIPS 24^m positions as priors to deblend 
the SPIRE maps. 

Although both the PACS and SPIRE catalogues have 
been produced in the same way, Magnelli et al. (2013) do 
not provide flux upper limits. In order to keep the priored 
FIR catalogues consistent with each other, we calculate up¬ 
per limits for the non-detections in the PACS catalogues of 
Magnelli et al. (2013) in a similar way to the upper limit 
calculation performed for the SPIRE priored catalogues of 
Swinbank et al. (2014). This was done by performing aper¬ 
ture photometry at thousands of random positions in the 
PACS residual maps and taking the 99.7th percentile of the 
distribution of the measured flux densities as the 3cr upper 
limit on the nondetections. To account for the fact that the 
deblending is more uncertain in regions of luminous sources, 
we calculated these 3cr upper limits as a function of the 
pixel values in the original maps (see Swinbank et al. 2014). 

2 The PACS catalogues for GOODS-N and GOODS-S are 
published in Magnelli et al. (2013). The catalogue for COS¬ 
MOS was created in the same way and is available online 
(http: / / www.mpe.mpg.de / ir/Research/PEP/DRl). 



Figure 1. X-ray (2-8 keV) luminosity (L2-8kev) versus redshift 
[z) for the X-ray sources in the GOODS-N, GOODS-S, and COS¬ 
MOS regions described in §2.1. Black centers indicate the X-ray 
sources without a direct hard-band detection (as described in 
§2.1). The vertical dashed lines indicate the 4 redshift ranges used 
in this study. The lower X-ray luminosity threshold (L2_8keV > 
10^2 erg/s) used to define our AGN sample is shown with the hor¬ 
izontal dashed line. The combination of the three fields enables 
us to explore the SFRs of AGN over three orders of magnitude in 
AGN luminosity. 


Consequently, this approach results in upper limits being 
higher for non-detected sources that lie near a bright source, 
when compared to non-detected sources in blank areas of the 
maps. 

Due to the fact that we are using MIPS 24/im pri¬ 
ored catalogues for the FIR photometry of our sources, 
any undetected at 24/im will not have FIR measurements 
in the published catalogues. Therefore for the 24/im unde¬ 
tected sources, we extracted the FIR photometry at the op¬ 
tical counterpart positions following the method described 
in Swinbank et al. (2014). Overall there are only 23 sources 
that are undetected at 24/im but have FIR counterparts, 
making up a very small fraction of our overall sample. 

In the MIR bands we also use the catalogues of 
Spitzer-IRAC 8/im observations as described in Wang et al. 
(2010), Damen et al. (2011), and Sanders et al. (2007), for 
GOODS-N, GOODS-S, and COSMOS, respectively, as well 
as Spitzer-IRS 16/im from Teplitz et al. (2011) for GOODS- 
N and GOODS-S. Since all the IRAC catalogues have their 
detections determined by the 3.6/im maps, and the 16/im 
catalogues have been produced with the use of 3.6/im pri¬ 
ors, they are all consistent with the deblended PACS and 
SPIRE catalogues described above. 

2.3 Redshifts and catalogue matching 

For our SED fitting analysis (see §3.1) we need matched 
catalogues containing X-ray fluxes, MIR-FIR photometric 
flux densities, and redshifts. To obtain the appropriate coun¬ 
terparts for each X-ray source, we matched the catalogues 
starting with the X-ray catalogues described in §2.1. We 
first match the positions of the optical counterparts of the 
X-ray sources to the MIPS 24/tm positions in the catalogues 
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Field 

AGN 

with spec- 2 : 

with 24/im 

GOODS-N 

177 

98 

137 

GOODS-S 

209 

128 

154 

GOSMOS 

1753 

914 

1151 

Total 

2139 

1140 

1442 


Table 1. Number of X-ray detected AGN in our parent sample 
(L2 —8keV > 10"^^ erg/s; z = 0.2-2.5) in each field, as well as the 
number of sources with a spectroscopic redshift and the number 
of sources with a 24 fim counterpart. 


of Magnelli et al. (2013).^ To choose the matching radii be¬ 
tween catalogues we measure the number of total matches 
as a function of radius and estimate the fraction of spuri¬ 
ous matches for each matching radius. The matching radius 
of the X-ray to the MIPS 24/tm catalogue for GOODS-N 
and GOODS-S was 0.8”, while for G-GOSMOS and XMM- 
GOSMOS it was 1”. This matching radius was chosen to 
maximise the number of matches while minimising the num¬ 
ber of spurious matches, with a ratio of spurious to true 
matches of 1%. Due to the way that the FIR catalogues 
were deblended, each MIPS 24/rm detected source also has a 
corresponding photometric measurement or flux upper limit 
for PAGS 70/im, lOO/rm, 160^m and SPIRE 250/rm, 350/rm, 
500/im (see §2.2). For the sources not matched to a MIPS 
24/rm counterpart we use the FIR data extracted at their 
optical counterpart positions, as described in §2.2. We then 
match to the IRAG, and to the IRS 16/rm catalogues for the 
two GOODS fields (see §2.2) using the same method. 

A necessity for this analysis are the redshifts of the 
X-ray sources. To allocate the redshift counterpart of the 
sources in GOODS-S and C-GOSMOS we make use of the 
spectroscopic and photometric redshift compilation by Xue 
et al. (2011) and Givano et al. (2012), respectively. We also 
added redshifts from Teplitz et al. (2011) for sources in 
GOODS-S when necessary. For the sources in GOODS-N 
we created our own compilation using catalogues of spec¬ 
troscopic redshifts from Barger, Gowie & Wang (2008) and 
Teplitz et al. (2011) and photometric redshifts from Wirth 
et al. (2004) and Pannella et al. (2009). Overall we obtained 
redshifts for 91.4% of the X-ray sources. 

In total there are 3297 X-ray sources covered by Chan¬ 
dra, XMM, and PEP/GOODS-i7 observations with a red¬ 
shift (see Figure 1). For this study we restrict this sample to 
redshifts of a = 0.2 - 2.5 and a luminosity range of L 2 _ 8 kev 
> 10^^ erg s“^, resulting in our parent sample of 2139 AGN. 
Of the parent sample 53.3% have spectroscopic redshifts and 
67.4% are detected at MIPS-24^m (see Table 1 for a sum¬ 
mary of the three fields). 


3 DATA ANALYSIS 

In this study we are interested in measuring the mean SFRs 
of galaxies, hosting an X-ray detected AGN, as a function 

® For the X-ray catalogues of CDF-S and C-COSMOS the optical 
counterparts are provided by Xue et al. (2010) and Elvis, Civano 
& Vignali et al. (2009). For the sources in CDF-N we use the 
catalogue of Barger, Cowie & Wang (2008). 


of the AGN luminosity and redshift. We use multi-band IR 
photometry, including photometric upper limits, to perform 
SED fitting for all 2139 X-ray detected AGN in our parent 
sample (see §2.3; Figure 1). For each source we decompose 
the contribution of AGN activity and star formation to the 
overall SED. This allows us to measure the IR luminosity 
due to star formation (Lir,sf), the key quantity for this 
study, which we can use as a proxy for SFR (e.g., Kenni- 
cutt 1998, Galzetti et al. 2007, Galzetti et al. 2010). In §3.1 
we outline the SED fitting procedure and describe the cal¬ 
culation of Lir,sf. In §3.2 we describe the method that we 
follow for the calculation of the average Lir,sf as a function 
of L 2 _ 8 kev (our tracer of the AGN luminosity) for the whole 
sample, where we include both direct Lir,sf measurements 
and upper limits. The calculation of these values thus al¬ 
lows us to investigate how SFR relates to AGN luminosity 
(Section 4). 

3.1 SED fitting procedure 

To calculate individual Lir^sf values for our sample we per¬ 
form SED fitting to the MIR and FIR photometry. In these 
bands there could be a contribution from both AGN and 
star formation, with emission from the AGN peaking at MIR 
wavelengths and dropping off at the FIR wavelengths (e.g., 
Netzer et al. 2007; Mullaney et al. 2011). Those factors make 
it important to decompose the contribution from both star 
formation and AGN to the overall SED so as to avoid an 
overestimation of the SFR measurement. In Figure 2 we give 
example SED fits to demonstrate our procedure. 

To fit and decompose the IR SED of our sources we 
develop the publicly available DecompIR code of Mullaney 
et al. (2011), and use the 8 - 500/rm data and upper limits 
described in §2.2. We use a set of empirical templates that 
consist of the mean AGN template and the five star form¬ 
ing galaxy templates originally defined in Mullaney et al. 
(2011), and extended by Del Moro et al. (2013) to cover the 
wide wavelength range of 3 - 10® /rm. We also include the 
Arp220 galaxy template from Silva et al. (1998) which serves 
as a sixth template to ensure that we are also covering the 
possibility of extremely dusty star forming systems. The ad¬ 
vantage of using a few, but representative, templates to fit 
the data is that we can avoid the degeneracy in the fitting 
procedure caused by a large number of templates. Further¬ 
more as many of our sources have limited photometric de¬ 
tections (with only one or two data points), it is sensible to 
keep the number of free parameters as small as possible. We 
note that the set of star forming galaxy templates described 
above covers a broad range of empirical shapes, including the 
large template library of Ghary & Elbaz (2001), as shown 
in Figure 2 of Del Moro et al. (2013), and the templates 
described by Kirkpatrick et al. (2012). 

In our fitting procedure the only free parameters of the 
fit are the normalisation of the star forming galaxy and AGN 
templates. Since there are two free parameters in the fit 
we require that the source has at least three photometric 
detections to simultaneously fit the AGN and star forming 
galaxy templates. When we have less than three photometric 
detections we can only derive upper limits on Lir^sf, as we 
cannot constrain the AGN contribution (see below). 

When a source is detected in three or more photometric 
bands we perform a series of fits following the method of Del 
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Figure 2. Examples of the four types of best fitting SED solutions, (a) A galaxy where the best fit (solid curve) is the combination of 
AGN (dashed) and star forming galaxy (dot-dashed curve) templates, (b) A galaxy where the best fit is that of a star forming galaxy 
template alone, with no AGN contribution, (c) A galaxy where the best fit solution is an AGN (solid curve) with no star formation 
contribution, in this case we calculated an upper limit on the star forming component shown by the grey dot-dashed curve, (d) A galaxy 
with only one photometric detection where we can only calculate an upper limit for the star forming galaxy templates, as shown here by 
the grey curve. In all four cases the blue data points are from Spitzer observations, while red data points are from Herschel observations. 
The filled circles are measured flux densities, while the empty circles with an arrow are the flux density upper limits. For each example 
we also give the SFR, X-ray luminosity, and redshift of the source. The wavelengths have been shifted to the rest frame, but the flux 
densities are in the observed frame. 


Moro et al. (2013). We fit the data in two steps: firstly we 
fit using each of the six star forming galaxy templates sepa¬ 
rately without including the AGN component, and secondly 
we fit again with each of the star forming galaxy templates 
in combination with the AGN template. We fit to the pho¬ 
tometric flux density detections, but use the available flux 
density upper limits to eliminate the fits which are above 
any of the upper limits. This procedure results in a maxi¬ 
mum of twelve models (the six star forming galaxy templates 
without an AGN and the six star forming galaxy templates 
with an AGN) to chose from. 

To determine the best fitting solution of the twelve pos¬ 
sibilities described above, we use the Bayesian Information 
Griteria (BIG; Schwarz 1978) which allows the objective 
comparison of different non-nested models with a fixed data 
set, and is defined as: 

BIG = -2 X InL -R k X InN (2) 

where L is the maximum likelihood, k is the number of free 
parameters, and N the number of data points. This method 
penalises against models with extra free parameters counter¬ 
balancing the fact that a model with more free parameters 
can fit the data better, irrespective of the relevance of the pa¬ 
rameters. This is an improvement over a simple Ax^ test or 
a maximum likelihood comparison that would tend to favour 
the model with more free parameters. For each source the 


BIG value is calculated for all of the different fits. The best 
fitting model will be the one which minimises the BIG value, 
its absolute value being irrelevant; however for one model to 
be significantly better than the others it needs to have a dif¬ 
ference in BIG value of ABIG 2. If ABIG ^ 2 then both 
models are considered equally valid (e.g. Liddle 2004). Our 
final best fit solution is the one with the lowest BIG value; 
however we only accept the AGN component as significant 
if the inclusion of it reduces the BIG value by ^ 2. In Figure 
2(a) we show a best fit SED that includes the AGN and star 
formation component, and in Figure 2(b) a best fit SED with 
only the star formation component. From the best fit SEDs 
we then measure the integrated 8 — 1000 fim IR luminosity 
of the star formation component (Lir,sf)- Furthermore, if 
multiple fits have BIG values equal to the minimum BIG 
value, we consider them equally valid and take the average 
of their derived Lir,sf. 

For sources detected in fewer than three photometric 
bands we can only calculate upper limits on Lir,sf, due to 
the insufficient degrees of freedom to calculate the AGN con¬ 
tribution to the IR luminosity. To calculate the upper limits 
of the normalisation of each star forming galaxy template we 
increase the normalisation of each template until it reaches 
one of the Scr upper limits, or exceeds the 3 ct uncertainty 
of a data point. We take the star forming galaxy template 
with the highest upper limit of Lir,sf as our conservative 
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upper limit for that source (e.g. see Figure 2(d)). Using the 
same method we also derive upper limits on the star for¬ 
mation contribution for sources where the best fit is fully 
dominated by the AGN (e.g. see Figure 2(c)). 

Due to the limited photometry and quality of the data, 
our procedure is not expected to significantly detect an AGN 
component in the IR SEDs of all sources. Indeed, the detec¬ 
tion of the AGN component in the MIR will be dependent 
on the relative ratio of Lir,sf over the IR luminosity due 
to the AGN (Lir^agn); for example, a source with a high 
ratio of Lir,sf over Lir_agn will not show strong evidence 
of an AGN component in its IR SED (e.g., see Appendix 
A of Del Moro et al. 2013). However, we note that if we 
force an AGN component to be present in the IR SEDs of 
each of our sources, our results of mean Lir^sf in bins of 
X-ray luminosity and redshift (see §3.2) only change within 
a ~5% level, which is smaller than the uncertainty of the 
mean Lir^sf results presented in §4.1. We also verified that 
our results were not sensitive to the choice of AGN template 
that we used, by refitting sources with two different AGN 
templates. One template is representative of low luminosity 
AGN, while the other template is representative of high lu¬ 
minosity AGN, as provided by Mullaney et al. (2011). The 
first template is “colder” than that used in our main anal¬ 
ysis, with less emission in the MIR and extended emission 
to the FIR wavelengths, and the second template is “hot¬ 
ter” , with most emission occurring at MIR wavelengths and 
a steep drop-off in the FIR (in agreement with the mean 
empirical templates of Quasars in the FIR; e.g. Netzer et al. 
2007). Between them, these two templates, encompass most 
dumpy-torus models (see Fig. 7 in Mullaney et al. 2011). In 
both cases our results of mean Lir^sf in bins of X-ray lu¬ 
minosity and redshift (see §3.2) only change within a ~10% 
level, which again is smaller than the uncertainty in the 
mean Lir^sf results presented in §4.1. 

Using our SED fitting approach we have a sample of 
2139 AGN with individual measurements (including upper 
limits) of Lir,sf. From our results for the whole sample there 
are 263 fits that required a significant AGN component in 
addition to star formation, 274 fits that required only the 
star forming galaxy template, and for 1602 sources only up¬ 
per limits on the star formation component could be derived 
due to limited photometry. 


3.2 Calculating average source properties 

For this study we aim to constrain the average star forma¬ 
tion rates of our X-ray AGN sample as a function of redshift 
and X-ray luminosity. A challenge for all studies using Her- 
schel FIR photometry is the low detection rate of individual 
sources (e.g., Mullaney et al. 2012b). In our sample we can 
only place upper limit constraints on the Lir^sf for many of 
our sources, i.e. 1612 out of the 2139 (75.4%) sources in our 
sample, due to the limited photometry or because they are 
AGN dominated. In our study we have placed conservative 
upper limits on the Lir,sf for the AGN for which it was not 
possible to directly identify the star formation component 
(see §3.1). In order to not bias our study to only the FIR 
bright sources we study the average properties of the whole 
X-ray selected AGN sample by using a Survival Analysis 
technique (e.g., Feigelson & Nelson 1985, Schmitt 1985) to 


calculate the mean IR luminosities with the inclusion of all 
of the upper limits (details below). 

We divide our sample in to four redshift ranges, a = 0.2 
- 0.5, 0.5 - 0.8, 0.8 - 1.5, and 1.5 - 2.5. For each redshift 
range we also divide the sample in to L 2 _ 8 kev bins deter¬ 
mined such that they included «40 sources in each bin. To 
ensure that all of the sources within the redshift range are in¬ 
cluded we allow the number to vary slightly, resulting in bins 
of 40 - 43 sources. For each L 2 - 8 keV “ 2 bin we calculate the 
mean IR luminosity due to star formation ((Lir^sf); see §3.1) 
and mean X-ray luminosity ((L 2 _ 8 kev); see §2.1). To calcu¬ 
late the (Lir^sf) values, with the inclusion of upper limits, 
we use the Kaplan-Meier product limit estimator (Kaplan 
& Meier 1958), a non-parametric maximum-likelihood-type 
estimator of the distribution function. We use the formula 
as described in Feigelson & Nelson (1985) for the estimation 
of the mean of a sample including the upper limit values. 
The advantage of this method is that it does not assume an 
underlying distribution. We will refer to this method as the 
K-M method for the rest of this paper. 

The main requirement for the use of the K-M method, is 
for the upper limit values to be randomly distributed among 
the measured values and independent of them. Due to the 
different types of upper limits that result from our fitting 
procedure (see §3.1) the upper limits on Lir,sf are indeed 
random enough for the use of this method. Furthermore, a 
K-S test on our L 2 - 8 keV - « bins, with a probability thresh¬ 
old of 1%, shows no evidence of the distributions of upper 
limits and measured values being drawn from different dis¬ 
tributions. This method also requires that the lowest Lir,sf 
value in each bin is a measurement and not an upper limit. 
For the 12 bins where this is not the case we follow the 
popular procedure amongst studies using this method, and 
assume that the lowest value is a measurement (e.g., Feigel¬ 
son & Nelson 1985, Zhong & Hess 2009). These 12 bins are 
randomly distributed with L 2 - 8 keV and redshift (see Table 
2) , and therefore do not affect our conclusions on the trends 
of (Lir,sf) with redshift, and L 2 _ 8 kev- 

Feigelson & Nelson (1985) use the K-M method to es¬ 
timate means with up to a censorship (i.e., the fraction of 
upper limits) of 90%, but argue that there can be a signif¬ 
icant bias in such cases. Additionally, a study by Zhong & 
Hess (2009) estimating the bias of this method for a wide 
range of distribution types, find that the estimated means 
are within a factor of 2 for up to 80-90% censorship levels. 
In our work we have imposed a limit of 90% censorship on 
our bins, and have discarded 7 bins with greater censorship. 
The median censorship level amongst the remaining 45 bins 
we have used in our analysis is ~73%, with 11 of them hav- 


Our SED fitting procedure provides upper limits for the cases 
where a source is: (a) MIR - FIR undetected; (b) MIR - FIR 
detected, but in less than three bands; (c) AGN dominated, i.e. 
the fit doesn’t require any contribution from the SF templates. 
In the case of (a) the upper limits are calculated by constraining 
the SF templates to the flux upper limits, while in the cases of 
(b) and (c) they are calculated by constraining the SF templates 
to the 3(7 flux errors or the flux upper limits. The fact that bright 
IR sources can meet the criteria of (b) and (c), in combination 
with the spatial dependence of the FIR flux upper limits, helps 
drive the similarity between the distributions of the Ljr sp upper 
limits and measurements. 
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ing a censorship of 80-90% (see Table 2). For the calculation 
of the uncertainty on the mean we use the bootstrap tech¬ 
nique, for which we take 10000 random samplings in each 
bin and recalculate the mean. We then take the 16th and 
84th percentiles of the overall distribution as the Icr errors. 
As discussed above, bins of high censorship levels could suf¬ 
fer from additional uncertainties of a factor of <2. However, 
when comparing to the results of the stacking procedure, we 
hnd that the two methods are consistent (see Appendix), 
and hence, we do not have concerns about the high censor¬ 
ship levels in our bins causing a significant systematic bias. 

We show our final results of (Lir^p) as a function of 
L 2 - 8 keV in Figure 3. In our plots, throughout Section 4, we 
also Include axes of SFR and AGN bolometric luminosity 
(Lagn) to help interpret the Lir,sf and L 2 - 8 keV measure¬ 
ments. We calculate Lagn from L 2 _ 8 kev by using the lu¬ 
minosity dependent relation of Stern (2015) to convert the 
L 2 - 8 keV to an AGN 6/rm luminosity density. We then mul¬ 
tiply this by a factor of 8 to convert the 6/im luminosity 
density to Lagn (following Richards et al. 2006). The SFRs 
were calculated from the (Lir^sf) with the use of the Kenni- 
cutt (1998) relation corrected to a Chabrier IMF (Chabrier 
2003). 


4 RESULTS & DISCUSSION 

In this section we present our results and explore the form 
of the relationship between the average SFR, (Lir^sf), and 
X-ray luminosity, L 2 - 8 keV 5 for our sample of 2139 X-ray de¬ 
tected AGN (see Section 2) . In §4.1 we present our results 
of average SFR (calculated from (Lir^sf)) as a function of 
X-ray (and bolometric) AGN luminosity for four redshift 
ranges within 2 = 0.2 - 2.5. In §4.2 we compare the SFR of 
the AGN to those of the overall star forming galaxy popula¬ 
tion, for a subsample of our sources with reliable host-galaxy 
masses. In §4.3 we compare our results to the predictions 
from two empirical models that connect AGN activity to 
star formation. 

4.1 Mean star formation rate as a function of 
X-ray luminosity 

The main focus of this paper is to determine the form of 
the relationship between the average SFR and AGN X-ray 
luminosity over 4 redshift ranges. The results of our anal¬ 
ysis as described in §3.2 are presented in Figure 3 and Ta¬ 
ble 2. In Figure 3 the data are colour coded by redshift 
where each point is the mean of «40 sources and the error 
bars correspond to the Icr of the bootstrap errors (see §3.2). 
We find that the (Lir_sf) (and hence (SFR)) increases with 
redshift, by a factor of ~3 between each redshift range, in 
agreement with both the observed evolution found for nor¬ 
mal star forming galaxies (e.g. Elbaz et al. 2011; Schreiber 
et al. 2015) and previous studies on AGN populations (e.g. 
Shao et al. 2010; Rovilos et al. 2012; Rosario et al. 2012; 
Mullaney et al. 2012b). However for the individual redshift 
ranges we hnd no strong correlation between (Lir^sf) and 
L 2 - 8 keV, a result inconsistent with that suggested by some 
other studies which have reported a rise or fall of (Lir,sf) 
at high X-ray luminosities (e.g.. Lutz et al. 2010; Page et al. 
2012; Rovilos et al. 2012; although see Harrison et al. 2012). 


We hnd that our results are in general agreement to 
those studies that stack the FIR data to derive SFRs using 
large numbers of sources (e.g., Harrison et al. 2012, Rosario 
et al. 2012); however our results have reduced scatter and 
reduced uncertainties on the AGN contribution to the IR lu¬ 
minosity. We look in more detail at how our results compare 
to those of stacking in the Appendix of this paper. Addition¬ 
ally, we compare our results directly to those of Rosario et al. 
(2012), who explore the average 60/rm luminosity (icLeo^m) 
values (as a tracer of SFR) in the same redshift ranges as 
our study, by stacking Herschel-PACS data. We use the av¬ 
erage diherence between icLeo^m and Lir^sf from our SED 
htting results, Lir,sf /icLeo/jm = 2.2, to convert the results 
of Rosario et al. (2012) to Lir_sf. In Figure 4 we plot our re¬ 
sults in comparison to those of Rosario et al. (2012) (hollow 
black symbols) and hnd broad agreement with our results 
both as a function of redshift and L 2 - 8 keV, although we have 
more L 2 - 8 keV bins and our results show less scatter. To com¬ 
pare to the highest L 2 _ 8 kev bins of Rosario et al. (2012) at 
the redshift ranges of z = 0.8 - 1.5 and « = 1.5 - 2.5 we cal¬ 
culate the (Lir^sf) for the hve highest L 2 _ 8 kev sources in our 
study in both of these redshift ranges (plotted in Figure 4 
with solid black symbols). We hnd that our highest L 2 - 8 keV 
sources are in agreement with those of Rosario et al. (2012); 
however, due to the very small number of sources in these 
bins (5 - 23 sources across both studies), we do not interpret 
them any further. 

To asses the contribution of the upper limits on the 
overall mean, we take an extreme scenario where all upper 
limits are assumed to correspond to zero values. We hnd that 
(Lir,sf) can drop by 0.2 dex (factor of 1 .6) at 0.2 < 2 < 0.5, 
by 0.3 dex (factor of 2) at 0.5 < 2 < 0.8 and 0.8 < 2 < 1.5, 
and by 0.4 dex (factor of 2.5) at the highest redshift range 
of 1.5 < 2 < 2.5. However, we note that the form of the 
observed flat relationship of (Lir^sf) with L 2 - 8 keV (Figure 
3) shows little to no change for all redshift ranges, in this 
extreme scenario. 

To test whether our results are consistent with a hat 
trend of (Lir_sf) with L 2 - 8 keV we show in Figure 4, as a hor¬ 
izontal grey line, the mean (Lir^sf) for each redshift range. 
Across all redshifts the data he within a factor of 2 of the 
mean. However, we hnd that the (Lir,sf) values of the most 
luminous AGN for all of the redshift ranges at 2 < 1.5 are 
systematically above the overall mean. To quantify the de¬ 
viation between the (Lir,sf) of the high L 2 - 8 keV bins to the 
rest of the data we make two hts; one to the two highest 
L 2 - 8 keV bins (with the exception of 2 = 0.2 - 0.5 where we 
use only the highest L 2 - 8 keV bin); and one to the rest of 
the luminosity bins in the same redshift range (see the grey 
dashed lines of Figure 4). We find an increase in (Lir,sf) 
by a factor of ~2 for the highest L 2 - 8 keV when compared to 
the lower L 2 - 8 keV bins in each of the redshift ranges with 
2 < 1.5. For 2 = 1.5 - 2.5 there is no signihcant difference 
in (Lir,sf) between the highest and lowest L 2 - 8 keV that we 
cover. We note that the systematic increase of (Lir^sf) at 
high L 2 - 8 keV values observed in the redshift ranges of 2 < 
1.5 does not correspond to a systematic increase of the red- 
shifts at high L 2 - 8 keV values (see Table 2). Thus the modest 
trends observed at the high L 2 - 8 keV are not driven by red¬ 
shift. We investigate the observed trends further in §4.3. 
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Figure 3. Mean IR luminosity due to star formation, {Ljpr gp), as a function of X-ray luminosity, {L2_8kev)) for four redshift ranges. 
Each L2-8keV bin has ~40 sources. We also give the corresponding SFR values using the Kennicutt (1998) relation corrected to a 
Chabrier IMF (Chabrier 2003), and the bolometric AGN luminosity Lagn calculated from L2_8keV using the luminosity dependent 
relation of Stern (2015). The errors on the (Lur gp) are calculated using the bootstrap analysis as described in §3.2 (see also §3.2 for a 
discussion on the additional uncertainties). 



Figure 4. {Lur^gp} as a function of (L2_8kev}i plotted in 
Figure 3 (also to be referred to for axis definitions). The hori¬ 
zontal grey lines indicate the overall mean Liu gp across all of 
the L2_8keV bins for each redshift range. The dashed grey lines 
indicate the mean Liu gp for (1) the one or two highest L2_8keV 
bins and (2) the lower L2_8keV bins for each redshift range (see 
§4.1). The black hollow symbols are the stacking results of Rosario 
et al. (2012), and the black filled symbols are bins of the high¬ 
est L 2 _ 8 keV sources from our study (we note that there are very 
few sources in these bins for both studies; see §4.1). Our results 
are broadly consistent with a flat relationship; however, for the 
redshift ranges with z < 1.5 the highest L 2 _ 8 keV bins are system¬ 
atically a factor of ^2 higher than the mean Liu gp. 


4.2 Comparing to the average SFRs of the overall 
star forming galaxy population 

Here we explore whether X-ray AGN have SFRs that are 
consistent with being selected from the overall star form¬ 
ing galaxy population. We compare the average SFRs of the 
AGN to the observed relationship between SFR, redshift, 
and stellar mass (M*) of normal star forming galaxies, which 
is defined as the “main sequence” of star forming galaxies 
(e.g., Noeske et al. 2007; Elbaz et al. 2011; Schreiber et al. 
2015; Speagle et al. 2014). To make this comparison we re¬ 
quire stellar masses for the AGN in our sample. We use the 
stellar masses from Ilbert et al. (2013) for the sources in the 
C-GOSMOS area. Since their analysis did not take into ac¬ 
count of a possible AGN component to the rest-frame UV 
to near-IR SEDs, we applied a colour cut to exclude sources 
for which there is likely to be significant AGN contamination 
to the SED. We only include AGN with rest frame colours 
U — V > 1 and V — J > 1 based on the analyses of Mul- 
laney et al. (2012a). This results in a subsample of primarily 
moderate luminosity AGN (L 2 - 8 kev^ 10"'"' erg s ') making 
up ~40% of the parent sample a.t z < 1.5, but only 26% of 
the parent sample at 2 = 1.5 - 2.5. For these sources, with 
a reliable M*, we calculate the (Liu.sf) as described in §3.2. 
Due to the reduced number of sources with masses we can no 
longer use bins of ~40 sources and we therefore reduce the 
number of sources required in each bin to 25. We show the 
(Lir.sf) as a function of L 2 - 8 keV, for the sub-sample with 
reliable M* values, in Figure 5. We note that this sub-sample 
have (Lir^sp) values consistent with the whole parent sam- 
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73 

3.2l°;®xl0-‘® 

86tl^ 

0 

1.86 

2.4x10'^® 

4.6xl0‘‘‘‘ 

85 

3.6lg;^xl0'‘® 

94tl^ 

1 

1.9 

3.0x10'^® 

6.3xl0‘‘‘‘ 

76 

4.2lg;^xl0-‘® 

11211^ 

0 

1.88 

3.8x10'^® 

8.5xl0-‘‘‘ 

81 

3.8lg;|xl0'‘® 

lOlt^l 

0 

2.02 

7.3x10'^® 

2.1xl0‘‘® 

83 

4.4lo;gXlO-‘® 

11611^ 

1 

1.94 

8.5x10'^® 

2.6x10"^® 

78 

3.7lo;|xlO'‘® 

98tl^ 

0 

1.95 

1.2xl0‘‘"‘ 

4.4x10"^® 

85 

3.8lg;?xl0-‘® 

loot"! 

0 

1.89 

1.4x10'^"' 

5.7x10-^® 

71 

5.6lp3XlO-‘® 

148tg 

0 

2.01 

2.1xl0‘‘"‘ 

l.OxlO'^® 

81 

3.2l°;®xl0-‘® 

86tl® 

0 

1.94 

2.6xl0‘‘"‘ 

1.4xl0‘‘® 

76 

5.7lg;^xl0'‘® 

15011® 

0 

1.91 

3.6xl0‘‘"‘ 

2.5x10"^® 

85 

3.llg;^xl0-‘® 

82II5 

0 

2.09 

8.0x10'^"' 

1.2xl0‘‘’’ 

83 

4.8lp2XlO‘‘® 

86lg 

1 


Table 2. The average redshift, X-ray luminosity, AGN bolometric luminosity, IR luminosity due to star formation, and SFR, for the data 
presented in Figure 3. The errors on the {Lur^sp} are calculated using the bootstrap analysis (see §3.2). We also provide the censorship 
level of each bin, and a flag indicating when the minimum value of the sources in that bin is an upper limit (when the flag has a value 
of 1), which can result to an extra uncertainty on the (Ljj^ gp) (see §3.2). 


pie (see Figure 5), with the exception of the a = 1.5 - 2.5 
range which appear to be systematically higher. 

We use the mean redshift and mean M* of each bin 
in Figure 5 to calculate the expected range in Lir,sf for 
“main sequence” galaxies using Equation 9 of Schreiber et al. 


(2015). The shaded regions, colour-coded by redshift, corre¬ 
spond to the range of Lir,sp covered by the main sequence 
galaxies at the mean redshift and mean M* of the sources 
in each bin; i.e. a scatter of 2 around the mean results from 
Schreiber et al. (2015). We also find that these results are the 
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Figure 5. (Ljr sf) a function of {L 2 _gkev) for the subsample 
of sources that have a reliable stellar mass (M*) measurement 
in Ilbert et al. (2013) (see §4.2; also see Figure 3 for the axis 
definitions). The grey solid lines are the means for each redshift 
range of the whole parent sample (see Figure 4). The shaded 
regions correspond to the expected range in Ljr sp for the overall 
star forming galaxy population at the mean redshift and mean M* 
of each bin as defined by Schreiber et al. (2015). For all redshift 
ranges the (Lijj^ gp) values of the AGN appear to be consistent 
with normal star forming galaxies. 


Figure 6 . The probability distribution of the Eddington ratio (A) 
for the three cases assumed in Figure 7 (i.e. two broken power law 
distributions with a faint end slope of a = —0.65 and a = —0.2, 
and a lognormal distribution with 0.4 dex dispersion; see §2.2 in 
Aird et al. 2013). This also serves as a schematic representation 
of the three distributions assumed for the Hickox et al. (2014) 
model, assuming that the shape of the distributions represent 
the variability function of individual AGN (see Section 4.3 and 
Section 2 of Hickox et al. 2014). 


same if we use the Elbaz et al. (2011) definition of the “main 
sequence”. We find that, for this sample of X-ray AGN with 
L 2 - 8 keV< 10"''* ergs“*, the (Lir,sp) in all redshift ranges 
with 2 : < 1.5 are consistent with that of star forming galax¬ 
ies of the same mean redshift and mass. This result agrees 
with the results of previous studies (e.g., Mullaney et al. 
2012b, Harrison et al. 2012, Bongiorno et al. 2012, Lanzuisi 
et al. 2015). However, for the redshift range of 2 = 1.5 - 
2.5 the (Lir^sf) is systematically at the higher end of the 
Lir,sf region covered by “main sequence” galaxies, which 
may be due, in part, to a bias due to the fact that only 26% 
of the parent sample at those redshifts have reliable masses, 
and these have systematically higher {Lir,sf) values than 
the parent population (see Figure 5). 

4.3 Comparing to empirical models 

As shown in Figure 4, the trend of (Lir^sf) ((SFR)) with 
L2-skeV (Lagn) is broadly consistent with being flat. This 
result may initially seem in disagreement with the results 
of studies such as Rafferty et al. (2011), Mullaney et al. 
(2012a), Chen et al. (2013), Delvecchio et al. (2014), and 
Rodighiero et al. (2015), which find a correlation between 
the average Lagn and SFR of star forming galaxies. How¬ 
ever, these studies start with a parent population of galax¬ 
ies for which they calculate the average Lagn, while in this 
study we start with a population of AGN for which we cal¬ 
culate the average SFR. It has been suggested that the vari¬ 
ability of AGN, taking place on smaller timescales to that 
of star formation, could flatten any intrinsic correlation be¬ 
tween the SFR and the Lagn when not averaging over the 
most variable quantity (i.e. by taking the average Lagn over 
bins of SFR; e.g., Hickox et al. 2014). To assess what could be 
the driver of the flat relationship that we observe, and if in¬ 


deed it is AGN variability that is driving its shape, we com¬ 
pare to two empirical “toy-models” that predict the (Lir^sf) 
as a function of Lagn- Firstly that of Hickox et al. (2014) 
and secondly, a model based on Aird et al. (2013) (also see 
Caplar, Lilly & Trakhtenbrot 2014 for a similar model). 

The empirical “toy-model” presented in Hickox et al. 
(2014) assumes that SFR is correlated to Lagn when av¬ 
eraged over timescales of the order of 100 Myr. To create 
the SFR distribution of the galaxies in their model, they 
assume the redshift dependent IR luminosity function (LF) 
from Gruppioni et al. (2013). In the model, the individual 
AGN are allowed to vary on short time scales on the ba¬ 
sis of an assumed Lagn /(Lagn) distribution, which serves 
as a tracer of the Eddington ratio distribution of individ¬ 
ual AGN in the absence of black hole masses (see details 
in Hickox et al. 2014). The fiducial model assumes that the 
distribution of Lagn /(Lagn) has the form of a Schechter 
function (broken power law form) with a faint end slope of 
a = —0.2 and a cutoff luminosity of Lent = 100 (Lagn) (see 
the dashed red curve in Figure 6 for a schematic of this dis¬ 
tribution). The model can then predict the average SFR as 
a function of instantaneous (i.e. observed) Lagn of a large 
population of simulated AGN. We ran the model for the four 
redshift ranges of this study and plot the results in Figure 
7(a) with solid tracks. The model successfully reproduces 
an increase of the (Lir,sf) with redshift, for a fixed range in 
Lagn, ^ and is in good agreement with the data at z = 0.2 
- 0.5; however, it fails to reproduce the trends observed for 

® The increase of the (Ljr sp) with redshift, for a fixed range in 
Lagn, could initially seem contradictory to the model’s original 
assumption of a correlation of SFR and the long term averaged 
Lagn- However, even though the increase of (Lir sf) will be 
accompanied by an increase in the long term averaged Lagn, 
there is not a significant difference in the range of instantaneous 
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(a) 



(b) 

Figure 7. (Ljr sp) a function of {L 2 _ 8 kev} four redshift 
ranges compared to model tracks from (a) Hickox et al. (2014) 
and (b) the extended Aird et al. (2013) model (see Figure 3 for 
the axis definitions). The solid lines in both cases show the pre¬ 
dictions of the models with their originally assumed Eddington 
ratio distribution. From the two models the one of Aird et al. 
(2013), which is based solely on observational data, is in better 
agreement to our results; however both models demonstrate how 
the flatness of the trends in our results are likely to be a con¬ 
sequence of the assumed Eddington ratio distribution, or AGN 
variability. We also investigate how different the trends are when 
assuming different Eddington ratio distributions in the two mod¬ 
els (i.e. two broken power-law distributions with a faint end slope 
of a = —0.65 and a. = —0.2 respectively, and a lognormal distri¬ 
bution for which we only show the tracks for 0.8< 2 <1.5 to avoid 
confusion; see Figure 6). The different assumed Eddington ratio 
distributions show significant differences in the predicted trends. 
See §4.3 for more details. 


between SFR and the long term (Lagn), or could be caused 
by the lack of an explicit Eddington limit in the model but 
rather a cut-off limit at high Lagn /SFR ratios (see Hickox 
et al. 2014 for details). We investigate how the predicted re¬ 
lationship varies with different variability prescriptions later 
^ on in this section. 

The second empirical “toy-model” that we have corn¬ 
s' pared to is based on Aird et al. (2013), which we extended 
tr. to make predictions for the relationship between AGN lu- 

"> minosity and star formation. This model uses the observed 
redshift dependent stellar mass function (SMF) of galax¬ 
ies (from Moustakas et al. 2013) in combination with the 
probability function of a galaxy of a given stellar mass and 
redshift hosting an AGN, based on measurements in Aird 
et al. (2012) for z <1. This model predicts the distribution 
of stellar masses, for which they correct to BH masses as¬ 
suming Mbh = 0.002 xM* based on Marconi & Hunt (2003), 
as a function of X-ray luminosity. In contrast to the Hickox 
et al. (2014) model they use an Eddington ratio distribution 
in the form of a broken power-law function with the faint 
end slope being steeper with a = —0.65 (see the blue curve 
in Figure 6 for a schematic of this distribution). ® In order to 
compare to our results we have extended the model to cover 
^ the same redshift range as that of our sample and convert the 
^ predictions of stellar mass to predictions of SFR. To achieve 
— this we adopt the measurements of the SMF by Ilbert et al. 
E (2013) at 2 ; = 1 - 2.5 as an extension of the Moustakas 
et al. (2013) SMF up to z = 1, and extrapolate the redshift- 
dependence of the probability of a galaxy hosting an AGN 
from Aird et al. (2012) to 2 > 1 (which is consistent with 
the 2 > 1 measurements from Bongiorno et al. 2012). Fur¬ 
thermore, we make the assumption that all of the AGN are 
hosted by normal star forming galaxies that lie on the “main 
sequence” as derived by Elbaz et al. (2011), which is moti¬ 
vated by the results of our study (see §4.2). ^ We convert 
from the model predicted stellar masses to SFRs, allowing 
for a scatter of 0.3 dex in SFR around the “main sequence” 
relation. In Figure 7(b) we present the resulting predictions 
of (Lir,sf) as a function of Lagn, plotted with solid lines, in 
comparison to our results for each of the four redshift ranges. 
The predicted trends of the mean SFR in this case are fiat 
for a wide range of Lagn, similar to our data, with a slight 
rise in (Lir,sf) at high Lagn (i.e., Lagn > ergs“^). 

On the basis of this modified Aird et al. (2013) model, the 
slight rise of (Lir^sf) observed in our data (see §4.1) may 
be driven by a small increase in the average masses of the 
galaxies hosting very luminous AGN. 

To first order, the data are better described by the 
extended Aird et al. (2013) model than the Hickox et al. 
(2014) model; see the solid tracks in Figure 7(b) compared 
to those in Figure 7(a). However, since the two models have 
assumed different Eddington ratio distributions (or, equiva- 


the higher redshift ranges. In particular, the normalisation 
of the predicted trends are too low compared to our data and 
the rise of (Lir,sf) with Lagn is much steeper than that ob¬ 
served. The steepness of the predicted (Lir^sf) trends at the 
highest Lagn could be a result of the enforced correlation 

LagNi across the simulated population, which is the quantity we 
effectively observe for an X-ray AGN sample. 


® We note that Aird et al. (2013) use an observed specific accre¬ 
tion rate distribution (i.e., Lagn /M*) which they convert to an 
Eddington ratio distribution. 

^ We note that there is evidence in optical studies of X-ray AGN, 
such as Azadi et al. (2015), that a small fraction of these AGN are 
hosted by non star forming galaxies; however, Azadi et al. (2015) 
find that these AGN appear to form a minority of the population 
and therefore we do not expect them to significantly affect our 
mean SFRs. 
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lently, Lagn /(Lagn) for the Hickox et al. 2014 model) we 
also explore how sensitive the results are to this assump¬ 
tion. We therefore also ran the models with a series of three 
different Eddington ratio distributions to understand how 
sensitive the predicted trends of (Lir^sp) with (Lagn) are 
on the assumed Eddington ratio distribution. We used (1) a 
broken power-law with a — —0.2 (i.e, the fiducial distribu¬ 
tion assumed by Hickox et al. 2014); (2) a broken power-law 
with a = —0.65 (i.e, the fiducial distribution assumed by 
Aird et al. 2013); and (3) a narrow lognormal distribution 
with a dispersion of ~0.4dex centred at an Eddington ratio 
of ~ 0.06, as defined by Kauffmann & Heckman (2009) for 
nearby AGN residing in star forming galaxies. These three 
distributions can be seen in Figure 6. 

In Figure 7(a)&(b) we show the three sets of tracks 
which correspond to the resulting trends of (Lir,sf) with 
(Lagn) for the different assumptions of the Eddington ra¬ 
tio distributions. A clear correlation between (Lir_sf) and 
(Lagn) is predicted for the lognormal distribution while, 
by comparison, the power-law models predict a much flat¬ 
ter relationship. With a change of power-law slope from 
a = —0.2 to a = —0.65, the normalisation of the model 
tracks increase and the trend becomes flatter. The different 
shapes of the model tracks are driven by the relative differ¬ 
ence between the low Eddington ratio slope and the slope of 
the low-mass end of the galaxy SMF (i.e., for M < M,, 
a ~ 0). The predicted correlation between (Lir,sf) and 
(Lagn) for the lognormal distribution is due to the narrow 
range of probable Eddington ratios. For the assumptions be¬ 
hind our models when assuming the lognormal distribution, 
most of the AGNs are accreting at a broadly similar Ed¬ 
dington ratio and therefore an increase in Lagn is predom¬ 
inantly due to an increase in stellar mass (and hence SFR 
since we assume the main sequence of star-forming galax¬ 
ies). By contrast, the steep low-Eddington ratio slope for 
the power-law models, when compared to the low-mass end 
slope of the galaxy SMF, allows for a broad range of Ed¬ 
dington ratios across a narrow range in stellar mass; i.e., 
there is a higher probability for an AGN of a given lumi¬ 
nosity to be hosted in a high-mass galaxy with a low Ed¬ 
dington ratio than a low-mass galaxy with a high Eddington 
ratio. Indeed, on the basis of the extended Aird et al. (2013) 
model, the population of low-to-moderate luminosity AGN 
(Lagn < 10^® ergs“^) predominantly reside in galaxies of 
similar stellar mass (M* ~ 10^®'®“^^M©), and thus similar 
SFRs, but with a wide range of possible Eddington ratios. 

Overall, our results suggest that the observed trends of 
(SFR) - Lagn are due to AGN being highly variable and re¬ 
siding, on average, in normal star forming galaxies. Similar 
results have also been found by hydrodynamical simulations 
that show that AGN variability can cause a flat trend be¬ 
tween Lagn and SFR (e.g., Gabor & Bournaud 2013; Volon- 
teri et al. 2015). The Eddington ratio distributions of AGNs 
are typically constructed to describe a population of AGN. 
However, as adopted in our models, they can also be under¬ 
stood as the distribution of Eddington ratios for an individ¬ 
ual AGN over time, and hence could be used as a variability 
prescription of the AGN (as originally adopted in Hickox 
et al. 2014). As is clear from Figure 7(a)&(b), the choice 
of Eddington ratio distribution plays a major role in the 
form of the predicted (SFR) - Lagn relationship. For ex¬ 
ample, our results are much better described with the use 


of a broken power-law Eddington ratio distribution with a 
faint end slope of a = —0.65, than with a narrow lognormal 
Eddington ratio distribution, which predicts a qualitatively 
different (SFR) - Lagn relationship to that found from our 
data. Thus, the (SFR) (or (Lir,sf)) - Lagn plane can be a 
useful diagnostic tool for placing constraints on the intrinsic 
Eddington ratio distribution of AGN (also see Veale, White 
& Conroy 2014). 


5 CONCLUSIONS 

We have created a large sample of X-ray detected AGN with 
FIR coverage and individual SFR measurements. Our sam¬ 
ple has a total of 2139 AGN at redshifts of z = 0.2-2.5, with 
10"^^ <L 2 - 8 keV< ergs“^. Using the available pho¬ 

tometry from 8-500/im we have performed individual SED 
fitting to all of the sources in our sample, and measure the 
IR luminosity due to star formation, Lir,sf. 

Our analysis has a number of key advantages over many 
previous studies: (a) the use of deblended source catalogues 
for the FIR photometry, which ensures better constraints on 
the flux density measurements and eliminates the overesti¬ 
mation due to blending and confusion of sources (see §2.2); 
(b) the use of photometric upper limits in the SED fitting 
analysis, which achieve better constraints on the fitted SEDs 
(see §3.1); (c) the decomposition of the AGN and star for¬ 
mation contributions to the FIR emission, which provides 
values of Lir,sf that are not contaminated by the AGN (see 
§3.1); (d) the calculation of upper limits on Lir,sf when the 
data are insufficient to identify the star forming component 
directly (i.e., not enough photometric data points, poor S/N 
data, or dominant AGN component), which allows us to es¬ 
timate the (Lir,sf) for all the sources in our sample avoiding 
the bias that could be caused by removing these sources (see 
§3.2). 

With the Lir,sf measurements for each source we de¬ 
rived the mean Lir,sf values ((Lir,sf); a proxy of the 
(SFR)) as a function of L 2 - 8 keV (a proxy of the AGN lumi¬ 
nosity; Lagn) in bins of ~40 sources, for the redshift ranges 
of 0.2 - 0.5, 0.5 - 0.8, 0.8 - 1.5, and 1.5 - 2.5. In comparison 
to previous studies, our results show less scatter amongst 
(Lir.sf) across the wide range of L 2 - 8 keV investigated in 
this study. Overall we found that: 

(i) The (SFR) increases by more than an order of mag¬ 
nitude from redshifts of 0.2 - 0.5 to 1.5 - 2.5, in agreement 
with previous studies on the redshift evolution of the SFR 
for the general star forming galaxy population. See §4.1. 

(ii) For each redshift range the (SFR) shows no strong 
dependence on AGN luminosity; however we note that for 
the redshift ranges of z ^5 1.5 the highest Lagn systems 
have (SFR) values that are systematically higher than those 
of lower Lagn systems by a factor of «2. See §4.1. 

(hi) For the ~40% of the sources within the GOSMOS 
area with reliable stellar masses, we compare their (SFR) 
to the “Main Sequence” of the overall star forming galaxy 
population. The X-ray AGN, at all redshift ranges, have 
(SFR) that are consistent with normal star forming galaxies 
at the same redshifts and masses. Due to a lack of secure 
masses for the high Lagn systems in our sample this result 
is restricted to moderate AGN luminosities (i.e., L 2 - 8 keV^ 
erg s-i or Lagn< 10^®'® erg s'^). See §4.2. 
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(iv) To qualitatively understand the flat relationship be¬ 
tween the (SFR) and Lagn we compared to two empirical 
“toy-models” that make predictions for this relation: Hickox 
et al. (2014) and an extended version of Aird et al. (2013). 
These models take mock galaxy populations and assign them 
with SFR values based on observed distributions, and in¬ 
stantaneous Lagn values based on an assumed Eddington 
ratio distribution. We And that the flat relationship seen 
in our data could be due to short timescale variations in 
the mass accretion rates, which, in combination with the 
relative shapes of the Eddington ratio distribution and the 
galaxy SMF, can wash out the long term relationship be¬ 
tween (SFR) and Lagn. See §4.3 

(v) We find that the predicted (SFR) - Lagn relationship 
is sensitive to the assumed Eddington ratio distribution. For 
example, both models predict a relatively flat relationship 
over all redshift ranges, assuming an Eddington ratio distri¬ 
bution of a broken power-law form with a faint end slope of 
a = —0.65, whilst with a log-normal distribution the pre¬ 
dicted trends are too steep to be consistent with our data. 
Therefore, the observed (SFR) - Lagn relationship appears 
to be a sensitive diagnostic of the intrinsic Eddington ratio 
distribution of AGN. See §4.3. 
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APPENDIX A: COMPARISON OF THE K-M 
METHOD TO THE STACKING ANALYSIS 
METHOD 

In this Appendix we compare our results using our SED 
fitting analysis and the K-M method that weused in this 
work (see Section 3), to those we would obtain with stacking 
analysis, a method commonly used in similar studies of star¬ 
forming and AGN galaxy samples. 

Following the method of Harrison et al. (2012) we 
stacked the SPIRE-250/rm maps at the X-ray positions of 
the sources of our sample in C-GOSMOS. We use the C- 
GOSMOS sample since it makes up most of our overall sam¬ 
ple and avoids issues that can arise when combining stacks 
of different fields with different depths. 

We bin the sample in bins of L 2 - 8 keV and reshift con¬ 
taining ~40 sources each, in the same way as described in 
§3.2 for the K-M method (in the redshift range of 2 = 0.2 
- 0.5 we use ~30 sources to allow for more than one bin). 
We show the stacking results in Figure Al, in comparison 
with the overall means of the K-M method results for each 
redshift range, as well as the results of Rosario et al. (2012). 
We find that our main results are consistent with the results 
we obtain when using the stacking analysis, and that both 
methods are in agreement with the results of Rosario et al. 
( 2012 ). 

This comparison demonstrates that our method for cal¬ 
culating the mean produces results consistent to the pop¬ 
ular method of stacking in the FIR. However, our method 
produces less scatter amongst bins, as well as smaller uncer¬ 
tainties on the mean values. This is likely due to the use of 



Figure Al. (Lir^sf) as a function of (L2-8kev} when stacking 
the SPIRE data at 250^m for the sources of C-COSMOS see 
Figure 3 for the axis definitions. We compare these results to the 
overall K-M means of our SED results (grey lines; see §4.1), and 
the results of Rosario et al. (2012). We find that our results are 
consistent to those obtained using the stacking analysis, however 
the K-M method’s results produce less scatter (see Figure 4). 


deblended FIR photometry, and the removal of AGN con¬ 
tamination, in our analysis, which are effects not taken into 
when stacking. 
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